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^ , Abstract. By using symmetry properties, the two-body Dirac equation in coordinate rep- 

' resentation is reduced to the coupled pair of radial second-order differential equations. Then 

the large-j expansion technique is used to solve a bound state problem. Linear-plus-Coulomb 
potentials of different spin structure are examined in order to describe the asymptotic de- 
04 I generacy and fine splitting of light meson spectra. 
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^ ■ 1 Introduction 

O 

^ I Two-body Dirac equations (2BDE), i.e., the Breit equations ^ and its generalizations [2113 

CLc m ini 13 El 111 ^1 ) are used frequently for the description of relativistic bound state problem, 

fH , especially in nuclear [111 I12j and hadronic |13| I14| I15| I16j physics. Apart from two free-particle 

Dirac terms, the 2BDE may include potentials which are local matrix-functions in the coordinate 
representation. This form provides an intuitive understanding of the interaction and may suggest 
a proper physical choice of the potential in phenomenological models. 

But the 2BDE are pathological if certain interaction terms are not treated perturbationally. 



B 

> 

^ ■ The set of radially reduced equations ^1 El El may possess non-physical energy-dependent 

poles at finite distance r between particles ^1 El . Correspondingly, an exact boundary- value 
problem becomes incorrect mathematically or improper for assumed physical treatment. 

Here we consider a possibility to avoid pathological peculiarities of 2BDE using a pseudo- 
perturbative technique similar to jl9l I2U1 OT] or 1/i expansions |22l I23j . These methods 
are applicable to the case of a strong coupling and are little affected by boundary peculiarities 
of the boundary-value problem. 

In our case natural expansion parameter is where j is the conserved total angular mo- 
mentum. After the radial reduction is performed, the 2BDE takes the form of the set of eight 
coupled first-order differential equations jSlIZI- Using a chain of transformations we reduce it to 
the pair of coupled second-order equations and apply the 1/j expansion technique. The method 
is applied to the potential model of meson based on the 2BDE. 

2 2-body Dirac equation and its radial reduction 

In the centre-of-mass reference frame the 2BDE has the form: 

{hi{p) + h2{-p) + U{r)-E}F{r) = 0, (1) 
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where F(r) is a 16-component wave function, 



haip) =CXa P + rUaPa = - i^a • V + TTlaPa, 



1,2, 



are Dirac Hamiltonians of free fermions of mass and U{r) is an interaction potential. If F{r) 
is presented in 4 x 4-matrix representation, the operators cxa and Pa act as follows: aiF = cxF, 
Q.2F = Fcx^ etc, where ex. and /? are Dirac matrices. 

The potential U{r) is the Hermitian matrix-function, it is invariant under rotation and space 
inversion transformations (so that the total Hamiltonian H = hi + h2 + U is too). Its general 
form is parametrized by 48 scalar function of r = \r\ jlUj . Of physical meaning are potentials 
admitting field-theoretical interpretation of interaction. In particular, potentials reflecting a spin 
structure of vector and scalar relativistic interactions are used frequently in potential quark 
models of mesons. We will consider such a model in the Section 5 using few examples of scalar 
and vector potentials known in a literature. In the present section the structure of potential is 
not essential. 

In order to apply a pseudoperturbative expansion method to the 2-body Dirac equation let 
us transform it to an appropriate form. 

First of all we perform a radial reduction. Following [21 [7j we put the wave eigenfunction 
F(r) of the total angular momentum j and the parity P into the 2x2 block-matrix form: 



F(r) = i 



isi(r)<^^(n) +iS2{r)4P{n) ti(r)(/)~(n) +t2{r)(p^{n) 
ui{r)(f)~ (n) + ti2(r)0+(n) ivi{r)(f)^{n) + iV2{r)<p^{n) 



(2) 



for the parity P = {—y^^ states, and into a similar form for the parity P = (— )-' states but 
with superscripts interchanged as follows: (^,0) <-> ( — ,-|-). Here n = r/r, the bispinor harmo- 
nics (/>^(n) corresponds to a singlet state with a total spin s = and an orbital momentum 
i = j, and (p^in), (j)~{n), </>"'' (n) correspond to triplet with s = 1 and i = j,j + l,j — 1. Then 
for i > the eigenstate problem (P) reduces to the set of eight first-order differential equations 
with the functions si(r), . . . , V2{r) and the energy E to be found. 

It is convenient to present this set in the following matrix form. Let us introduce the 8- 
dimensional vector-function: X(r) = {si(r), S2(r), ti(r), . . . , V2{r)}. Then the set of radial equa- 
tions reads: 



X(r) = EX{r) 



(3) 



where the 8x8 real matrices H(j) and \l{r,E,j) = G{j) / r + m-\-\}{r, j) — E\ possesses properties 
H"'" = — H, V"'" = V, the diagonal matrix m = diag(m+J, m„J, — m_J, — m+J) (here J is 2x2 unit 
matrix and m± = mi it 7712) and j- and P-dependent matrices H(j), Q{j) are constant (i.e., 
free of r), and matrix-potential U(r, j) comes from interacting term of the equation Q. For the 
case j = components S2 = i2 = ^2 = ^2 = so that the dimension of the problem @ reduces 
from 8 to 4. 

It turns out that rankH = 4 (2 for j = 0). In other words, only four equations of the 
set (jnj are differential while remaining ones are algebraic. They can be split by means of some 
orthogonal (i.e., of 0(8) group) transformation. In new basis we have 



X(r) 



Xi 
X2 



H 



J(2) 




where J^^^ is the symplectic 4x4 matrix. Thus we arrive at the set 



2j(2)x;+VnXi + Vi2X2 = 0, 



(4) 
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V21X1 + V22X2 = 0. (5) 

Eliminating X2 from @ by means of © we get a differential set for the 4-vector Xi 

+ V^(r,^, j)| Xi(r) = 0, where = (Vn - \l^2M22^2i) /2 

while X2 then follows from the algebraic relation X2 = — \/^2^V2iXi. 

The elimination of X2 causes non-physical energy-dependent singular points (apart of r = 
and physical singularities of potentials) in matrix elements of V^. 

Now we present the 4-vector Xi in 2 -|- 2 block form, 



Xi(r) 



$1 
$2 



Vn V12 
V21 V22 



eliminate then $2 and arrive at the second-order differential equations for 2- vector <I>i: 

!i{E)^i = { (|: + ^12) [^22]"' - V21) + Vnj = 0. (6) 

The matrix V22 is diagonal for all potentials considered in Section 5 (and many other ones). 
In these cases we can perform the transformation: 



1.1 = $i/VV22, ^ 



V22i:vv 



22 



providing for the operator L the form which is as close as possible to 2-term form: 



Z{r,E,j), 



d 
dr 



(7) 



here W(r,E,j) is a symmetric 2x2 matrix, is 2 X 2 symplectic matrix and {•,•}+ denotes 
the anticommutator. 

We are going to apply the 1/j expansion method to the equation (O. In many physically 
interesting cases the function Z vanishes or it is negligible at j large. Thus the wave equation has 
a 2 X 2 matrix 2-term form which is convenient for application of the method. In other cases the 
third term of the operator ((T)) contains a first-order derivative via off-diagonal matrix elements 
only. This form is tractable too, but with more tedious calculations. We do not consider such 
equations in this paper. Before proceeding further, we study a simpler example of a single 2-term 
relativistic equation. 

3 Todorov equation via 1/i method 

Here we consider the Todorov-type equation describing the relativistic system of two interacting 
scalar particles in the centre-of-mass reference frame [241 1251 126j : 

{p2 + U{r, E) - b{E)} ^(r) = 0. 

Here p = — i V, the quasipotential U (r, E) depends on energy E of the system, and the binding 
parameter b{E) is the following function of E, 



h{E) = \e'^ -\{ml + mfj + \{ml-ml) /E"^, so that E{h) = ^ ^/ml + b. (8) 



a=l 
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The corresponding radial equation takes the form 



W{r,E,i)'^^>{r)=Q, (9) 

where ^ is the angular momentum quantum number, and 

W{r, E, £) = U{r, E) + £{£ + I) /r^ - b{E) . (10) 

Let us consider motion of the system in the neighbourhood of classical stable circular orbit. 
Given ^ > 0, the radius = rc{i) of the stable circular orbit and the corresponding energy 
Ec = Ec{i) satisfy conditions: 

W{r„E„£)=0, dW{r„E„i)/dr, = (11) 

and d^W{rc,Ec,£)/dr^ > 0; here dW{rc, Ec,£)/drc = dW {r, E,£)/dr\ r^rc etc. 

\E—Ec 

One puts r = Tc + Ar and E = E^ + AE where Ar and AE are small in some meaning, 
and expand the function W{rc + Ar,Ec + AE,i) in power series with respect Ar and AE. 
Then due to the conditions ()11() the leading terms of this expansion represent the harmonic 
oscillator potential and other ones are anharmonic terms. If the conditions (|11() hold for any 
large value of i it is possible by renormalization of Ar and AE to single out in the equation Q 
the ^-independent harmonic oscillator problem and anharmonic perturbations which disappear if 
i ^ oo. This is the idea of 1/^ expansion method. Application of pseudoperturbative techniques 
of this type ^1 [201 1^ |22l 122] to our case meets two peculiarities: the equation Q represents 
a nonlinear spectral problem, and an exact solution of the equations Hll|) may appear to be 
unknown or too cumbersome for practical use. Thus we modify the technique. 

Let us introduce the parameter A = l/^/I which is small at £ large. Since the exact form of 
the functions rc{i) and Ec{i) is unknown in general, we first determine asymptotics rc ~ roo(A), 
be = b{Ec) ~ 600(A) at A — > which may be found much easier. Then the functions rc(^) 
and Ec{i) can be presented in the form: 

rc(A) = roo(A)p(A), b,{X) = b^{X)fi{X), 

p(X) = 1 + Ap(^) + AV^'^ + • • • , /i(A) = 1 + A^i(i) + aV^'^ + • • • , (12) 

where expansion coefficients = 1,2,... (and thus the analytical functions p(A) and 

/i(A)) can be found, step by step, from the conditions: 

W{p,fi,X) = 0, dW{p,fi,X)/dp = (13) 

and d'^W{p, p, X)/dp^ > 0; here the dimensionless function W{p,p,X) is constructed by the 
direct use of (\12\i in 1)10(1 and normalizing in order that W{p, p, A) to be regular at A — > 0, 

W{p, p, X) = AV^T^ [rooP, E{b^p), l/A^] . 

Now we go to the dimensionless variable r — > and spectral parameter b{E) e, 

r = roo(A)[p(A) + Ae], b = booiX) [fi{X) + Xh] , (14) 

in terms of which the equation Q takes the form 

{^-^^(^'^'^)}V'(O = (15) 
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with 

v^(e) = ^[roo(p + Ae)], 

and 

w{^,e,X) =W{p + X^, fi + Xh, X). (16) 

If the functions p{X) and fi{X) satisfy the conditions (fT!^ . the equation (P3|) is nonsingular at 
A — > 0. This is true even if we use the first-order approximate solution to (|13p in H14|) . 

p{X) = l + Xp^^\ /i(A) = l + A^W. (17) 
Indeed, using the notation dW^'^ydp = lim dW /dp. = dW/dp{p = l,p = 1,A = 0) etc. we 

A — ^0 

have 

^«.({,e, A) = -^W [p(\) + A{, m(A) + X\ A] 

Singular terms are absent if the following set of equations holds: 

iy(o) = 0, dW^^'^/dp = 0, (19) 

Besides, zero-order terms which are linear in ^ disappear if 

Notice that the equations 1)19^ and (|2U() - H2ip represent the conditions (|13() in the zeroth and 
first orders of A, respectively. Thus the equations (|19j) hold identically and (|2flj) - (|2H) are linear 
equations with p(^) and //^^^ to be found. 

In zero-order approximation the equation (|15|) reduces to the harmonic oscillator problem 



^ip^ + f^e-u- cu^f } -0(0 = (22) 
with 
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The higher-order terms in the expansion (|18() can be considered as perturbations of the oscillator 
problem (|22jl . They depend, in general, on the spectral parameter e and can be taken into 
account by means of the perturbative procedure is appropriate to this case. Otherwise the 
treatment is similar to QUI 1^ 1^ OH] . 

The eigenvalues in zero-order approximation e„,, = [u){2nr + 1) + y]/ where n,. = 0, 1, . . . 
is a radial quantum number, are to be corrected by means of higher orders of perturbative 
procedure. Then, using of the 2nd equation of (|14() in ((HJ gives us the energy spectrum. 



4 Breit-type equation via 1/j method 

At this point we return to the radial 2BDE in the form L{E)^i = 0, where the 2x2 matrix 
operator ^{E) is given by equation ^ with the last term neglected. Let us put 



$1 



^-1 

^2 



where ^'i and ^'2 are components of $1. Then the equation can be presented as a pair of 
coupled ordinary second-order differential equations: 

^^i{r)-Wi{r,E,j)-^^{r)=Y{r,E,j)m2{r), (26) 

-^^2(r) - W2{r,E,j)^2{r) = Y{r,E,j)^i{r). (27) 

We will treat this system perturbationally using the pseudosmall parameter A = l/-y/J- 

Let us suppose for a moment that the right-hand side of the system (|26 ]) -(|27 |) can be ignored, 
so that these equations decouple. Then we can apply to each of the equations the scheme of the 
Section 3. We define radii and energies of circular orbits by means of the conditions: 

W,i,M)=,. i?!l%^>0, i = l,2. 

or or-^ 

Then we single out asymptotics of these functions of A by means of the relations: 

ri{X) = 7'joo(A)pj(A), hi{X) = 6ioo(A)^fi(A), 

p,{\) = 1 + \pf^ + AVf ^ + • • • , /i(A) = 1 + Xpf^ + AVf ^ + • • • , 
and, using the relations 

r = ri^{\)[pi{X) + h = hi^{\)[p,{\) + i=l,2 (28) 

we reformulate the equation H26() in terms of the dimensionless variable .^1 and the spectral 
parameters ei while the equation (|27|) - in terms of .^2 and €2- Finally, we perform expansion of 
the equations into powers of A and solve them separately. 

Now we are going to take actual coupling of the equations ()26|) and (|27() into account. First 
of all, we note that the variables ^1 and ^2 are not of one another, and the spectral parameters ei 
and 62 are also not independent. Thus we should choose common variables in both equations. 

Let us first choose C = Cii ^ = ^i- Then the set (pH|) - (|27j) reduces to the form: 

< (0 - e, x)MO = e, X)MO, (29) 

i^UO - ^MC, e, x)MO = e, x)MO, (30) 
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where 

V'ilO = *ihoo(pi + Ae)], i = l,2, (31) 

w,{i, e, A) = X^rl^Wi [n^{pi + A^), i?(6ioo(/ii + A^e)) , l/A^] , (32) 

y{i, e, A) = X^l^Y [n^{pi + XC),E{h^{fi^ + A^e)) , l/A^] . (33) 

The functions H32() -()33 p are regular at A — > 0. Moreover, the general structure of the function 
wi is the same as that of w in the Section 3 (see equations (|18j) ). In particular, wi = O(A^). 
Thus the equation (|29() is similar to (|15|) (but with non-zero right-hand side). It admits similar 
expansion in A. 

On the contrary, the function W2 may have a different behaviour at A ^ 0. Here we consider 
three cases. 

1. Let r2oo 7^ ^^loo and 6200 / ^loo- Then W2 = 0(A~"), n > (except perhaps very special 
examples which we do not consider). In this case one can solve formally the equation (|3()j) in 
favour of ip2iO as follows: 

^2 = - f 1 - ' = - E f-^)" (34) 

V / W2 :^q\'"J'2 0^ J W2 

This representation leads to the loss of solutions for ip2 which are not analytical in A and thus 
have nothing to do with the perturbation procedure. The use of (|34j) in the r.h.s. of (|29|) permits 
us to eliminate "02 from H29|) and thus to obtain a close wave equation for ipi . The structure and 
treatment of this equation are the same as those of the equation p5|) . Moreover, it is obvious 
from (|34|) that at least the zero- and the first-order terms of 1^2 vanish. Thus the r.h.s. of (|29|) 
does not contribute in lower orders of perturbation procedure. In zero-order approximation we 
have the oscillator problem. 

2. Let r2oo = ?'ioo and 6200 = ^loo but p2 - Pi = 0(A) and /i2 - Pi = 0(A). Then W2 = 0(A). 
Since \^ /w2 = 0(A) the perturbative treatment H34() of the equation H30() is still valid. The only 
difference from the case 1 is that the r.h.s. of equation (|29|) may contribute in the first order 
of A. 

In both the above cases we used the dimensionless variable .^1 and obtained a closed eigenstate 
equation (which we will reference to as the problem 1) for the wave function V'i(?i) and the 
spectral parameter ei. We can proceed with the variable ^2 and obtain the problem 2 for the 
function 1^2(^2) and the parameter €2- One might be inclined to think that both problems 1 
and 2 are equivalent and lead to the same spectrum (in terms of energy E). Actually, different 
problems complement one another. This is evident from equation (|28j) leading to the relation: 

1 f^lcXD 1 , f^ioo ,1 

^^-^^=A^i^^^-^^ri^-'r^- 

Indeed, in both 1 and 2 cases |e2 — ei| — > 00 if A ^ 0. It does mean that an arbitrary energy 
level E calculated by means of eigenvalue e^''^ of zero-order oscillator problem 1 (with the use of 
equations (|28jl and ^) cannot be obtained by means of any finite eigenvalue €^2"^ the problem 2 
and vice versa. Higher-order corrections to e^''^ (or e^2^ ) are small and do not change qualitatively 
this picture. Thus different problems generate different branches of the energy spectrum of the 
original set of equation. In this respect the following special case differs essentially from the 
previous ones. 

3. Let r2oo = ^^loo and 6200 = ^loo but P2 - Pi = 0(A") and P2 - Pi = 0(A"), n > 2. Then 
W2 = O(A^). Both equations (|29jl and (|3U|) have similar structure and should be treated on the 
same footing. Use of common variables 4, e defined by (|TH) and H17() is appropriate to this case. 
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In the zero-order approximation we obtain the coupled pair of wave equations (on the contrary 
to the cases 1 and 2 where we had a single wave equation). In physically meaningful cases (of 
Section 5, for example) they have the form: 

{<f/de + Ate - - Lo^e] MO = XMO, (35) 
{cf/df + Ke-V2- ^2^2} V^2(e) = XMO, (36) 

where x = li™ V = const, and parameters Mj, k and lv are related to the functions Wi (i = 1, 2) 

A— >0 

by the equation of the type of (|23|) . (|24|) and (|25|) . The equations (|35jl . (|36|) can be evidently 
reduced to the pair of similar equations but with parameters vi = {z^i+i/2±Y^ {^i — V2)^ + 4x^}/2 
(i = 1,2) and x = 0. Thus they become split equations of the form l|22j) . The eigenvalues e 
corresponding to the first and second equations are separated by finite constant v\ — v^- Thus 
the corresponding states mix in higher orders of perturbation procedure. 



5 Application: Regge trajectories of mesons 

Here we apply the pseudoperturbative treatment of 2BDE in meson spectroscopy. 

It is known |27j that spectra of heavy mesons are described well by the nonrelativistic potential 
model with QCD-motivated funnel potential ii(r) = u;(r) + uc(r), where 

uc{r) = -a/r, a = 0.27, (37) 
ui{r) = ar, a = 0.25 ^ O.SGeV^ (38) 

The Coulomb part (|37|) of this potential describes a nonrelativistic limit of the vector one- 
gluon exchange interaction while the linear part ()38() is suggested by the area law in the lattice 
approximation of QCD and has presumably scalar or scalar-vector nature. 

Description of light meson spectroscopy needs application of appropriate relativistic models. 
Most of them are related to the string theory. From the theoretical viewpoint the most interesting 
are QCD-motivated relativistic models embracing properties of both heavy and light mesons. 
Such models should reflect the scalar-vector structure of interaction and should lead to funnel- 
type potential in the nonrelativistic limit. 

A natural candidate for the relativistic potential model is the2BDE with a short-range vector 
potential and a long-range scalar one. At least three general structures of vector potential are 
used in the literature, 

C/v(r-) = nv(r), (39) 
C/v(r) = {l-ai-a2}nv(r), (40) 
U^,{r) = {l - \oLi ■ Q!2}nv(r) + \{n ■ ai){n ■ a2)ru'^{r), (41) 

with Uv{r) = uc{r) or another short-range potential; here u'{r) = du{r)/dr. The potential (|39j) 
is only a static part of vector interaction (see |Sj). The relativistic vector field kinematics is 
taken into account in the potential (|40jl (see |16| |H]) which, for the Coulomb case, was first 
proposed by Eddington and Gaunt UHl • In the generalization H41|) of the Breit potential P 
retardation terms have been added |Hj. Two different scalar potentials, 

Us{r) = PiP2Us{r), (42) 
C/s(r) = i(/3i+/?2)ns(r), (43) 



come from different couplings of scalar mediating field with fermionic fields. The first one (|42|) 
arises from the Yukawa interaction (see [H]) while the second one (jlHjl corresponds to so called 
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"minimal" coupling The latter and also two following potentials can be treated as static 
approximation of various QFT-motivated scalar quasipotentials |3U1 1311 0^ I17j : 

C/s(r) = l(l + /?i/32K(r), (44) 
C/s(r) = i(l + /3i)(l + /32K(r). (45) 

The perturbative treatment of Breit-type equations has been used for calculating a fine 
splitting in spectra of heavy mesons jl41ll6j . Light meson spectra are essentially relativistic and 
need a nonperturbative statement of the problem which is inconsistent because of non-physical 
singularities of radial equations. To avoid these difficulties in numerical calculations one is 
forced to invent sophisticated potentials and impose rather artificial boundary conditions |15j . 

Using the pseudoperturbative treatment of 2BDE with different combinations of potentials 
H37|) - ()45() we obtain analytical expressions for meson mass spectra and estimate a role of general 
structure and input parameters of potentials in the model. We consider mass spectra of lightest 
mesons (containing u and d quarks only) and try to reproduce their following general features: 

i) Mass spectra of light mesons fall into the family of straight lines in the (i?^, j)-plane known 
as Regge trajectories. 

ii) Regge trajectories are parallel; slope parameter a is an universal quantity, a = 1.15GeV^ = 
(4^4.5)a. 

iii) Nonrelativistic classification of light mesons as (n^*'''^£j) -states of quark-antiquark system 
is adequate; i.e., radial quantum number rir = n — i — 1 enumerates leading (n^ = 0) and 
daughter (n^ = 1,2,...) Regge trajectories, spin s = 0, 1 corresponds to mass singlets and 
triples etc. 

iv) Spectrum is £s-degenerated, i.e., masses are distinguished by i (not by j) and n^. 

v) States of different i possess an accidental degeneracy which fact causes a tower structure 
of the spectrum. 

vi) Hyperfine ss-splitting is relatively small, about 5 6 % of cr. 

For this purpose we use the nonrelativistic potential function (|37j) and (|38jl in vector and 
scalar potentials of different spin structure and calculate pseudoperturbative spectrum 

in zero-order approximation. Classification of states then is done using singlet-triplet properties 
of large-large component of wave function @ in the nonrelativistic limit. 

If the vector short-range interaction is ignored and scalar potentials ()42() - H45|) are used with 
Us{r) = ar the pseudoperturbative mass (i.e., energy) of meson in zero-order approximation has 
the following form: 

E\ = ka[l+\ + r]{nr + ^)] + Cm+V2^ + Jim^ - 62mim2 + 0{l/VI), (46) 

here m+ = mi+m2, and /c, 77, 5i, 52, >i are dimensionless constants depending on the potential 
chosen. 

Four families of energy levels Ei {i = ^,0,—,+) form trajectories in the (i?^, £)-plane which 
are nearly straight. Indeed, C = 2 for all the potentials considered and rest masses nia 
(a = 1,2) of lightest mesons are small compared to the energy scale ^/a. Thus the parameter 
^m_|_-y/a determining a curvature of trajectories is small. Parameters 5i and 82 determine a com- 
mon shift of all the trajectories and are not important for the present discussion. Below we 
discuss the calculated values of parameters /c, x and 77 determining the slope of trajectories and 
their degeneracy properties. 
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In the ()42() case A; = 4 so that the slope a = ka matches quite well to that of property ii); 
■q = 2 causes accidental degeneracy typical for the harmonic oscillator; but x = 4 leads to 
j-dependence of energy (not ^-dependence) so that the /s-degeneracy is absent. 

In the case k = A and r/ = 2, so that the slope and the accidental degeneracy are the 
same as in the ()42() case; x = 4 — 3\/2 ~ —0.243 provides an approximate Zs-degeneracy, with 
accuracy 6%; the splitting is of order of the actual ss-splitting (see property vi)). 

In the ()44() case k = x = 3^/3 K, 5.196, r] = ^/2,K, 1.732; none of these values match well to 
properties ii)-vi). 

In the (113) case k = \/23 - \/T7(7 + \/T7)2/128 w 4.2 provides the best fit of a to that 
of property ii); rj = (VT? — 3) \/ 102 + 26\/T7 /8 ~ 2.03 leads to nearly precise oscillator-like 
degeneracy, with accuracy 1.5%; x = provides exact Zs-degeneracy. 

Taking into account the vector short-range interaction (one of potentials H39() - ()41|) with 
tiv(r) = uc{r)) results in a parallel shift of Regge trajectories. The value of the shift is of 
the order aa, it depends on the vector potential chosen and is different (in general) for different 
trajectories Ei (i = ^, 0, — , -|-). 

It has been proved in the framework of single-particle Dirac equation a possibility of confi- 
nement by means of vector and equally mixed vector-scalar long-range interactions I33^ 1341 . 
We examined these cases in 2BDE approach using different vector potentials (|39j) - (|4H) with 
liv(r) = ar. Corresponding zero-order pseudo-perturbative spectra have a form similar to ()46() . 
The difference is that k = ki = 8-^ 12 is two times or more larger than desired, and is different 
ior i = A,0,—,+ (i.e., trajectories are not parallel). 

6 Summary 

The Breit equation and its generalizations (2BDE) possess non-physical singularities. In some 
cases these points lay far from the physically important domain but they make a boundary 
problem incorrect or physically improper |17[ I16 j . In order to avoid this difficulty and to use 
the 2BDE in the relativistic bound state problem, especially for the case of strong coupling, we 
develop the 1/j expansion method. 

The method is based on the large- or large-£ techniques applicable to the radial Schrodinger 
equation. In our case the 2BDE is reduced to the coupled pair of quasipotential-type equations 
which structure causes principal modification of known techniques. Other changes are related 
to the fact that the equations represent a nonlinear spectral problem with cumbersome quasipo- 
tentials. 

We apply this pseudoperturbative method to the 2BDE with the linear -|-Coulomb potential of 
different scalar-vector structure. In all cases in the zero-order approximation it was obtained the 
Regge trajectories which are linear asymptotically. Linear potentials of two scalar structures H43|) 
and (|45|) which was discarded in ^\ as nonphysical (because of singularities in 2BDE) reproduce 
well in our case general properties of light meson spectra. In particular, the slope a = ka oi 
light meson trajectories fit well to the experimental value if the parameter a is taken from the 
nonrelativistic potential model [^Tj- The third linear potential of ^3 (with no singularities in 
2BDE) does not match to experimental data. 
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